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We address a study of glueball masses in the confining regime of ST/ (2) in D = 3 using an algo- 
rithm inspired by the multi-level scheme. Our method, which exploits the locality of the action to 
achieve high precision results, is based on a technique already used for compact QED, and gen- 
eralises it to the non-Abelian case. We discuss the main features of this method, in comparison 
with other algorithms that have been used in similar studies. 
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1. Introduction 

The basic idea underlying lattice calculations of the glueball spectrum in pure gauge theory is 
very simple: given a zero-momentum lattice operator with quantum numbers J PC , the connected 
correlator at large time separations At is expected to behave like: 

(<!>Hto + At)<j>(to)) coim * |(vac|0|O)| 2 e- EoAt , (1.1) 

where E$ is the energy of the lightest state |0) that couples to (p. The operator must have a suf- 
ficiently large overlap with the wave-functional associated to the required channel, and must be 
sufficiently smooth over the physical length scales characterising the wave-functional itself; this 
can be achieved by means of standard variational methods [1,2] and through "smearing" [3] and 
"blocking" [4] techniques, which allow to obtain precise results for SU(N) gauge theories [5-8]. 
However, these calculations are technically non-trivial, due to the fact that the connected corre- 
lators are typically very small quantities, affected by large relative errors at the time-like separa- 
tions where one expects to extract a signal for the mass; a solution was proposed in [9, 10], using 
anisotropic lattices with a fine spacing in the time-like direction, in combination with an improved 
action. 

In this work we discuss another method which is expected to improve the result precision 
in measurements of glueball correlation functions: following the idea proposed in [11], we use 
a multi-level algorithm to calculate the observable expectation values from products of sublattice 
averages; in particular, for the scalar channel, we consider a method which involves (discretised) 
time-like derivatives of the disconnected correlator. Since this method has been successfully used 
for compact QED [12], we are encouraged to test it in a non-Abelian model, focusing on SU(2) 
lattice gauge theory in 2 + 1 dimensions. Another study of the Yang-Mills spectrum using the 
multi-level algorithm can be found in [13, 14]. 

2. The LW multi-level algorithm . . . 

Liischer and Weisz [11] proposed a multi-level algorithm which achieves exponential error 
reduction in measurements of Polyakov loop expectation values in quenched LGT; given the two- 
point Polyakov loop correlation function: 

G(r) = (P*(r)P(0)) , with:P(x)=tr(jlU (t,xfj , (2.1) 

P*(r)P(0) can be factorized in a product of two-link operators T(t): 

P*(r)P(0) = |l]Tw} , ^ith:T(t) a p Y5 = Uo(t,ry ap U (t,d) r8 . (2.2) 
^ t ' aayy 

whose expectation values can be estimated to a high degree of precision from sublattice averages 
— conventionally denoted with the [. . .] notation — in time-like slabs of given thickness: the LW 
algorithm thus amounts to a variant of the multi-hit method [15], applied to pair of links; the scheme 
can be generalised to several nested sublevels. Due to locality of the gauge action, the different 
sublattices are independent of each other, provided the space-like link variables at their boundaries 
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are kept fixed; the thickness of the sublattices must be such that they are still in the confined phase. 
The value of the P*(r)P(0) product in the full lattice is then obtained multiplying the sublattice 
averages of the T's, and finally the two-point Polyakov loop correlator can be evaluated averaging 
over a sample of full-lattice configurations: in practice, after every (standard) update of the whole 
lattice, the algorithm performs a cycle of updates in the nested sublattices, keeping their boundaries 
fixed, to measure the [T(t)T(t + a) . . .] sublattice averages. 

The LW algorithm increases the result precision by damping the uncertainty due to short- 
distance fluctuations, through the sublevel averages: in a theory with non-vanishing mass gap, the 
algorithm leads to exponential error reduction with respect to standard simulations; if (at the lowest 
level) the whole lattice is subdivided in n s time-like slabs, then from n measurements in each slab 
(which require an amount of computer-time comparable with n full-lattice measurements in the 
standard approach) one effectively achieves a precision corresponding to n' h measurements in the 
whole lattice. Clearly, the efficiency of the algorithm depends on some parameters which can be 
tuned by optimisation. 



. adapted to study the glueball spectrum 



To generalise the LW algorithm idea to the study of glueballs, the correlator in eq. ( |1 ■ 1| ) has 
to be expressed in a factorized form, in terms of quantities that can be measured to a high degree 
of precision from sublattice averages; the most obvious choice is to consider the zero-momentum 
component of <p, evaluated in two time-slices belonging to different slabs. However, for the ++ 
channel one has then to subtract the non-vanishing VEV of the operators, keeping into account 
only the measurements that have been used in the two-point function, as it was remarked in [14]. 
Here, we follow a slightly different strategy, first used in the case of compact U(l) LGT in four 
space-time dimensions [12], considering time-derivatives of the full correlator: 

4-4-^Hh)^(to)}~-am 2 (e- n < t, - to) +e-'< N '- n+to) ) . (3.1) 
at\ oto \ / 

On the lattice, the expression above can be approximated in terms of finite differences, to read: 

F(h -t ) = V£v£(0l"(fi)0(lb)) ~ 2a +<? -»(M-ri+'D) J (i _ coshm ) . (3.2) 

Then we can use the LW algorithm to measure: 

v{v{(^(/iMlb)) = <[0 f (fi + 1) - <S>Hh)] [0(f„ + 1) - 0(fo)]) • (3-3) 

Notice that the quantity in eq. ( p\3| ) is expressed as a product of only two factors, whereas in 
the original formulation of the LW algorithm [11] the Polyakov line correlator was completely 
factorized in a (generally) larger number of factors. 



4. Algorithm details 

We consider a cubic lattice with sizes N x = N y = N t = 48, and use the standard plaquette action: 

S = ^(\-^trU n ) , with:/3 = ^2 ; (4.1) 
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Method 


arriQ++ 


y 2 

A- red 


Standard algorithm 


0.802(58) 


1.29 


Multi-level algorithm 


0.767(31) 


0.57 



Table 1: Comparison of the estimates for the lowest-lying mass in the ++ channel, obtained from runs 
requiring a similar amount of CPU time. 



the configurations are updated using a combination of the Kennedy-Pendleton heat-bath [16] fol- 
lowed by three overrelaxation [17] steps. 

The results discussed here have been obtained from simulations at j3 = 9.0, which corresponds 
to a lattice spacing a ~ 0.074 fm [18]; in the measurements, we focused on the J PC = ++ channel, 
considering a family of operators obtained from rectangles with sizes ranging from 1 to 4, built 
from "smeared" variables [3]. 

In the practical implementation of the multi-level algorithm, one needs to tune some parame- 
ters according to optimisation criteria (see [11]), including the number of nested average sublevels, 
the thickness of the slabs, and the number of measurements in the sublattices. As it was pointed out 
in [19], it is quite remarkable that the LW algorithm does not necessarily require a large number 
of nested levels: for most applications, a two-level algorithm proves to be highly efficient (and our 
experience in the present study confirms this empirical observation); in any case, more complicated 
schemes, with many nested levels, can be implemented numerically in a simple, versatile and ef- 
ficient way using recursive functions. An important parameter in the original LW algorithm is the 
thickness of the sublattices: for best efficiency, it has to be set to the minimum value compatible 
with the fact that the sublattice system is still in the confined phase, and a reasonable lower bound 
is half the inverse deconfinement temperature at the given value of j8 [18,20]. Finally, a "thumb- 
rule" was suggested in [13], estimating the optimum number of sublattice measurements in a slab 
of width L slab to be of the order of exp (moL s i ab ), where mo is the mass of the lightest state compatible 
with the quantum numbers of the operator which is measured. 

In the preliminary results of the study that we present here, we kept the mentioned criteria in 
mind; a systematic study of the parameter choice is among the subjects of work in progress, and 
we refer the reader to a forthcoming, complete paper [21] for a more exhaustive discussion. 



5. Preliminary results, observations, and comments 

The information emerging from the preliminary results of this study shows that the multi- 
level algorithm improves the estimate of the lowest-lying glueball mass: we observed that the 



estimate for the lowest mass in the ++ channel obtained from a fit to eq. (3^) is affected by a 
smaller uncertainty, with respect to a standard algorithm requiring a similar computational effort 
(see table [I]); also, the multi-level algorithm allows to explore a larger range of time separations 
(see figure |l|), where the signal obtained in standard simulations is already drowned in statistical 
noise. 

We also observed that our choice of the parameters in the study of the ++ channel has proved 



effective; in particular, the systematic error induced by the finite differences in eq. (3.2) did not 
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Figure 1: Results of the multi-level algorithm for the discretised double derivative of the correlator in the 
++ channel F(At) versus the time-slice separation in lattice units; the fitted mass reads: ma — 0.767(31). 

bias the estimate of the lowest glueball mass with respect to the results available in literature [5] — 
at least within the (limited) precision of the preliminary study presented here. 

Finally, we would like to point out that another advantage of using the multi-level algorithm 
in the study of the glueball spectrum (where, typically, one measures a large number of operators) 
is the reduced amount of data storage which is required, with respect to to standard simulations; 
this aspect, which was also remarked in [14], is simply due to the fact that the correlation matrices 
are built using the sublattice averages of the operators in a given channel, rather than the single 
measurements. 

6. Conclusions, perspectives and acknowledgements 

The analysis in progress is confirming that the LW multi-level algorithm can be successfully 
applied to numerical investigations of the glueball spectrum in non-Abelian lattice gauge theories: 
a conclusion in agreement with the results of similar studies. 

In particular, we would like to highlight the following aspects: 

• efficient error reduction with respect to standard simulations; 

• reduced data storage; 
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• versatility: the algorithm can be used for a generic gauge theory, and it can be combined with 
other methods which improve the quality of the results. 

A natural generalisation of our algorithm would be to use anisotropic lattices with a fine spac- 
ing along the time-like direction [9, 10]: this would enhance the data resolution at intermediate 
distances, with the further bonus of reducing the systematic error induced by the discretised time- 
derivatives used for the ++ case. 

Complete results of this work will be published in [21]. 

The author thanks Philippe de Forcrand, Ferdinando Gliozzi, Miho Koma, Yoshiaki Koma, 
Biagio Lucini and Antonio Rago for useful discussions, and acknowledges support received from 
Enterprise Ireland under the Basic Research Programme. 
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